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Abstract. Within the real-time approach [121 [E], the current across a quantum dot 
which is tunnel coupled to two leads at different chemical potentials is calculated in 
two steps: One determines the reduced density matrix of the quantum dot by using 
the "density matrix kernel" , and then the current by applying the "current kernel" to 
the reduced density matrix. If one multiplies the tunneling Hamiltonian by a coupling 
parameter "id", then everything, including the kernels and also the current, becomes 
a function of w. In the time space, the kernels have the clear structure of a power 
series in w, so one can speak about "orders" of the kernels. All odd orders and the 
constant terms vanish. Intuitively, one would guess that one can take the kernels of 
2n-th order in order to obtain the current of 2n-th order, i.e., the coefficients of the 
Taylor expansion of the current in w = up to the order 2n. I want to show that this 
is true. I can address only the single impurity Anderson model (SIAM) [3]. 
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1. Introduction 

The SIAM quantum dot [5] is a quantum dot with only four possible states: It can 
be empty ("0"), occupied by an electron with spin o ( "cr") or be in the state "2", 
i.e., filled with two electrons of opposite spin. If the quantum dot is tunnel-coupled to 
contacts at different chemical potentials, then a current can flow between them. Within 
the diagrammatic real-time approach [T2l |8] it is given by the equation 



where K c is called the "current kernel" and is the stationary reduced density matrix 
of the dot. The latter contains information about the stationary probabilities of finding 
the dot in the possible states. The current kernel is originally obtained as a function 
of time. The object of Eq. ([I]) is the Laplace transform of this function of time. One 
applies the map 



Jo 

to K c (t), where A is a positive number, and takes the limit A — > 0. 

The stationary reduced density matrix is determined by solving the quantum master 
equation in the stationary limit: 



The operator H Q is the Hamiltonian of the isolated dot. It is diagonal in the introduced 
states of the quantum dot: H Q (a) = E a a, where E a is the eigenenergy of the state 
a = 0,cr, 2. The object K{\ = 0) is the "density matrix kernel"; its structure is very 
similar to that of the current kernel and, initially, it is a function of the time; application 
of the Laplace transform to K(t) and the limit A — > yield K(X = 0). 

A possible approach to the problem of determining the current is perturbation 
theory: One can introduce a coupling parameter which expresses the strength of the 
tunneling coupling and consider the kernels and finally also the current as a function of 
this parameter. Then one calculates the Taylor series of u I(w)" in w = up to an order 
z as high as possible. The resulting polynomial I^ z '(w) of degree z is a function about 
which one can hope that it is a good approximation for I(w) in the case of small values 
of w. Since the current is calculated via the kernels, it is natural to try to obtain the 
Taylor series of the current by calculation of the Taylor expansions of the kernels. 

The relation to Ref. [7] is that in this work approximations for the kernels of a 
different kind were calculated: They are not perturbative since they are not Taylor 
expansions. The contributions to the kernels are given by what we call diagrams. All 
diagrams within a special class of diagrams which we refer to as the "DSO" were summed 
up. This class contains diagrams of all orders. The diagram selection has also been 
discussed in Ref. [10]. Moreover, e.g. in Ref. [8], another diagram selection called the 
resonant tunneling approximation (RTA) has been used. 



<I l > oc =Tr{K c (X = 0)p oo } 



(1) 




(2) 



= -[ Poo ,H Q ]+K(\ = 0)p, 



(3) 
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2. Hamiltonian 



The Hamilton operator of the SIAM quantum dot can be written as 



H Q = U d\.d^ d^di + E d\d a ^ 



(4) 



where d\ (d a ) is the creation (annihilation) operator of the one-electron level with spin 
a. The operator acts on the four-dimensional complex vector space spanned by the 
subsets of the set of the two one-electron levels, {t)4}- I abbreviate them by "0", "a" 
and "2". The eigenvalues are E = 0, E a = E, E 2 = IE + U. 

The contribution of the contacts to the total Hamiltonian is assumed to be given 



This is the usual choice of the Hamiltonian of the leads. 

Finally, there is the tunneling Hamiltonian, which expresses the possibility that 
electrons may tunnel from the leads to the quantum dot or vice versa [I] . I assume here 
the conventional tunneling Hamiltonian: 



where c[ kcr (Qko-) is the creation- (annihilation) operator of the electron level in the lead 
/ with wave vector k and spin a. For simplicity, I will assume that the coefficients TJk CT 
of the tunneling Hamiltonian are independent of spin, XJko- = T^. 

3. Taylor expansion of the kernels 

3.1. The kernels in the time space 

The kernels are obtained by taking the sum of all diagrams, the diagrammatic series is 
infinite. A possible way of writing down the density matrix kernel in the Laplace space 
was explained in Ref. [6]. It is obtained from the conventional diagrammatic language 
(e.g. Ref. [12]) by grouping the diagrams according to their topology. Also in the time 
space, the density matrix kernel can be written down in this way. For any operator x 
acting on the vector space spanned by the quantum dot states one obtains: 



by 





(hermitian conjugate), 



(5) 




2n 





n=l k=0 





CP'?) J: consistent with (p,q) 
(-l) {n+k) sign{p,q) |/3 >< a 2n \ < (3 k \x\a k > 
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J dk 1 Z h ...J 



dk n Z, 



i--Pi<qi<k 
i:fc+l<Pi<gi 

II /J^W*) < /VJ^JA, >< "a-J^V* > 

i:pi<k,k+l<qi 




n+T2+...+T2 ra _i<| 



2n-l 

n 



exp <^ -itj [ A£y + v Pi £ k*m ~ Yl v Pi £ h^i 

The meaning of the appearing objects is defined in Ref. [6], but I briefly explain 
the terminology here and also included Fig. [TJ 2n is the order of the diagram, n 
the number of its tunneling lines; k is the number of quantum dot states appearing 
on the left-hand side of the diagram; the (3i, an are arbitrary quantum dot states; I 
used the bra-ket notation. Z i; <7j is the lead- and the spin index of the i-th tunneling 
line; (p, q) is a set of n pairs of numbers between one and 2n and "is" the set 
of the tunneling lines; always Pi < qi. J is an increasing sequence of "intervals" 
J(l) C J(2) . . . C J{2n — 1) C {0, 1, ... , 2n} and defines the time ordering of the 
diagram; J(j) := J(j) \ minJ(j). The "pair formations" (p, q) originate from Wick's 
theorem [3] and are required to be "consistent " with the sequence J of intervals in such 
a way that the diagram is irreducible: For every j G {1, . . . , 2n — 1} there must be at 
least one i G {l,...,n} such that Pi G J(j),qi 4- J{j) or y i ce versa. This condition 
becomes graphically manifest in the property of Fig. [T] that for any horizontal square 
bracket representing one of the intervals J(j) there is at least one tunneling line which 
leaves the region enclosed by the square bracket. 

kj is the wave vector of the electron level attached to the tunneling line i, Z\ i 
the number of allowed wave vecors per volume in the first Brillouin zone of lead /j. 
The integrations go over these. By "/j " the Fermi-Dirac distribution/ one minus the 
Fermi-Dirac distribution of lead I is denoted. The operators Acrk : = T^d a contain the 
tunneling coupling; the integral with respect to the time is performed over the set 

t 

v ri, . . . , r 2n _i) : n > 0, ri + . . . + r 2n _i < - 

n 

which has the size ^n-i)! > ^ ne v p* are s iS ns which one can assign to the chosen sequence 
of quantum dot states, depending on their particle numbers. The Tj are the lengths of 
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a 



a, 



■4 






^ x T 2 k 2 o 2 




< 





a, 



corresponds to: 




Po Pi-^-[P 2 / a 2 ]-^-a 3 -J-a 4 



Figure 1. An example of a fourth-order diagram, written in the conventional 
(above) and an alternative (below) form. 2n = 4 in this diagram, moreover, 
k = 2; the pairs (pi,qi) appear in the lower version of the diagram as pairs 
of horizontal lines, i.e., pairs of numbers between 1 and 4 once these are 
counted from the left to the right. The pairs in the example are (1,4), (2,3). 
The sequence of intervals is given by J(l) = {2,3}, J(2) = {2,3, 4}, J(3) = 
{1,2,3,4}. Indeed, the pair formation and the sequence of intervals are 
consistent in the sense that the diagram is irreducible. 



time intervals between subsequent tunneling events. One multiplies them with the 
energy differences of the states with respect to which the total density matrix is off- 
diagonal during these intervals; is the energy of the one-electron level in lead I and 
wave vector k: AEj = E g . r , .. — E a T , ., . 

Because the SIAM has only the quantum dot states 0,7,},, and 2, the expressions 
containing the matrix elements of the annihilation operators of the quantum dot can, 
as far as they are not zero, always be written as 

Next, the integration with respect to the wave vector can be replaced by an integration 
with respect to the electron energy, as sketched in Ref. [Ij. One obtains "coupling 
functions" ai(e) which are, roughly speaking, proportional to the square of the tunneling 
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coefficients and to the density of electron levels in lead I around e. Upon introducing 
atf{s) '■= fi ± \e)ai{e), one can rewrite the part of the diagram beginnning with the 
integration over the Brillouin zones as 



± 



J de\ - J de n 




n ^ } ^) n n 

i:pi<qi<k i:k+l<pi<qi i:pi<k,k+l<qi 

2n-l 

, n 

Tl+T2 + ...+T2n-l< ^ j = l 

exp I -it 3 I AEj + ~ v P> £i f ( 6 ) 

In the exponent, the sums with respect to j and % can be swapped; the integration 
with respect to the times can be swapped with the energy-integrals, one obtains then 
Fourier transforms. To realize this plan I subdivide the set of the tunneling lines into 
four subsets (of {l,...,n}): 

h ■= {t:p l >k + l}, 

I2 ■= |i : 3j such that pi G J{j), q% ^ J(j) and: pi < k, qi > k + 1 j , 

h ■= {i ■ qi<k} , 

I4 := |i : 3j such that pi J(j), qi G J(j) and: pi < k, qi > k + l| . 

{1, ...,n} is the disjoint union of h, I2, h, and I4. The expression ^ reads in 
terms of Fourier transforms: 

2n-l 



± 




e X p(-i^r,A^) (7) 



n^R-*') k E 
n^(^' , )(-" W e 

™GJ 3 \ j.p-^J(j),q t &J(j) 



n^K"" 1 ) U, e 

i:ie/ 4 \ r.pigj(j),qi£J(j) 
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3.2. Expansion of the transform K(X = 0) 

In the time space, the kernels have the diagrammatic expansion 

oo 
71=1 

where w is the coupling parameter and is the contribution of the diagrams of order 
2n. The task is to prove that the expansion survives the transformation to the Laplace 
space ([2]), even in the limit A — > 0. To this end, I must show that for "sufficiently small 
values of if", i.e., for w < w , where w > 0, the series of integrals 



E 



00 

(2n) , 



n=l Jo 

is finite. 

Hence, back to the time integral, expression ([7]): The integral with respect to 
the r-s can be calculated by fixing a time t' between zero and the upper bound, t/h, 
performing the integration over the set where the sum of the Tj-s equals this time, 
and then integrating with respect to r'. In this way, the integrated kernel is related 
to its time derivative. Moreover, the appearing Fourier transforms are bounded by an 
exponential decay: |J r (. ..)(r)| < a(w) exp(— c\r\) = w 2 a w= i exp(— c\t\); this is the case 
at least if one makes a simple assumption about the coupling function, for example 
Lorentzian shape. In a subsection below I will show why this is the case. Making use 
of this inequality and with the definition 

Nj := number of elements of |i : p, L G J(j), q% J (J) or vice versa J 

one obtains the estimate 



w 



.In 



2» 



KS:\( t )\<i/h w ^j2 E 

k=0 /3o,...,/3fc,«fe,...02 

E E 

(P'i) J: consistent with (p,q) 
1 1 



E 



l\<j\,...,l n a r , 
n 



h y/2n - 1 



Tl + ...+T2r, 



exp -cJ2 T 3 N J 

i=t/» V j J 



The notation a w= i =: a was used. So far, I chose first an arbitrary pair formation (p, q) 
and then summed over all time orderings J which are consistent with the pair formation 
in the sense of irreducibility. These two things can be swapped. Taking the integral of 
the right hand side over the positive time axis, one obtains: 
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w 



*|Afflw|<i/» ^ E E E 

fc=0 /3 ,...,/3 k ,a k ,...a2„ l 1 a 1 ,...,l n cr n 

E E 

J (p,q): irreducible in J 

2n-l 

n ^ w 



3=1 3 



3.3. Combinatorial intermezzo 
Now, an estimate for the summand 



2n-l 1 

e n ^ < 9 ) 



(p,q): irreducible in J 



3=1 3 



is needed. The objects which define the quantity are an arbitrary ordering J of 
the numbers 1,2, ...2n, and a pair formation which is irreducible in J. For any 
j G {1, . . . , 2n — 1}, Nj is the number of pairs which are separated by the "interval" J (J), 
i.e., Nj is the number of is in {1, .. . ,n} with the property that Pi G J(j),qi 4- J(j)> 
or pi ^ J{j),Qi G Without loss of generality I can assume that the ordering J is 

just given by the natural way of counting. The sum ^ does not depend on J. One can 
view Nj = N(p, q)j as a map 

N(p,q): {l,...,2n-l}^{l,2,...} 

with the property that, with every step, the value of the function changes by ±1 and 
that N(p, q){2n — 1) = N(p, q)(l) = 1. The sum (|9| can be rewritten as 

function(N(p,q)) = function(N). 
(p,q): irreducible N (p,q)-N( P ,q)=N 

With this I express that there are in general different irreducible pair formations 
(p,q), {p',q') with N(p,q) = N(p',q'). One can take the sum over all irreducible pair 
formations by first fixing a map N, summing over all (p, q) with N(p, q) = N, and finally 
taking the sum with respect to N. All information about iV is contained in the subset 

S + :={l}u{j:N(j)>N{j-l)}. 

There are 2 2n subsets of {l,...,2n}, and, as a consequence, less than 2 2n ways of 
choosing the map N. 

Now, I want to fix the map N and think about 

2n-l 



2 function(N) = J] — - ^ 

(p,q):N(p,q)=N j=l U J (p,g):7V(p,g)=7V 

I define S~ := {1, . . . , 2n} \ S + . For sure, 
{l,...,2n} = S + US~. 
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I imagine that the numbers 1, . . . , 2n are listed in the natural order from the left to the 
right. If the pair formation (p, q) has the property N(p, q) = N, then the partner of 
a number / in S + can always be found to the right of /, the partner of a number /' in 
S~ can always be found to the left of I'. Thus, the number of possibilities to choose a 
partner for I' G S~ is always bounded by the number of those elements of S + which 
can be found to the left of I refer to this number as Z + {1'). Going through the 
elements of S~ from the left to the right, one notices that, for the first element /i, one 
has really Z + (l\) possibilities to choose a partner, for the second element, however, only 
— 1, for the third Z + (l 3 ) — 2, and so on. One finds that the number of different 
pair formations with N(p, q) = N is 

n n—1 

n^ + ^)-(^-i))=n( iv (y+ 1 )- 

i=\ i=l 

Finally, I used the characterization N(j) = \S + D {1, . . . ,j}\ — \S~ fl {1, . . . ,j}\ of 
the map N. One obtains: 

2n-l n-1 

funcUon(N) = J] — J] (N(l t ) + 1) 

(p,q):N(p,q)=N 3=1 ^' i=l 

< ff N{k) + 1 < 2-'. 
I can summarize the result of the combinatorial intermezzo by the estimate 

2n-l 1 

(p,q): irreducible in J i~ l 



3.4- Back to the expansion of the transform K(X = 0) 

This intermediate result can now be inserted into the inequality (|8|. For a fixed value 
of k, the number of different time orderings J is bounded by 2 2n_1 , since, everytime 
the interval is enlarged, there are at most two possibilities for this (towards the left 
or, alternatively, towards the right). The sums with respect to the spins and the lead 
indices as well as the sum with respect to the quantum dot states can be estimated 
in an elementary way, just by the number of possible choices. Finally, the number of 
possibilities to choose k is 2n + 1. One arrives at an inequality of the form 



POO 

^0 



< const! const n w 



2n 



with const, const' > 0, independent of n. One can conclude that K(\ = 0, w) is given 
by a power series in w 2 with positive radius of convergence. The coefficients are indeed 
K^ 2n \\ = 0). The treatment of the current kernel is analogous. The combinatorial 
intermezzo was necessary, since a crude estimate by the number of irreducible pair 
formations would have induced the presence of factorials. This causes problems, while 
the emergence of powers of some constant does not. 
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3.5. Why do the Fourier transforms decay exponentially? 
I assume functions 

a t {x) = — ^, (10) 

I _l_ I X-XQ 



CO 

1 



X — X± 

1 + e c i 

of a real variable x. I want to show that the Fourier transform of the product of 
the two functions, "a^~", decays exponentially. (The Fourier transform of on decays 
exponentially; it follows then that the Fourier transform of ai(l — fi), too, decays 
exponentially.) Without loss of generality I can assume that X\ — 0, c\ — 1. 

Basically, my idea was that the Fourier transform of a product of two functions can 
be written as the convolution of the Fourier transforms of the two functions. In order 
to realize this, I defined for A e]0, 1[ the function fx by 

fx(x) := ./;(•'' )< ' • 

For a fixed value of A, f\ decays exponentially. Moreover, there is the pointwise 
convergence 

f x a t -> a+ (A ->■ 0), 

and the convergence also has an integrable upper bound, since ai is integrable. I 
conclude: 

(Faf) (y) = lim F(fxai)(y) 

for every single value of y. 

Because of the relation (J 7 /)^/) = \f ) (y)> which can be shown via 

integration by parts [5], one knows that the Fourier transforms of fx, of on, as well as of 
their product, decay quickly, for example: quadratically. (I used here the smoothness of 
ai.) This ensures that the three Fourier transforms are integrable. Because the Fourier 
transform is invertible, the equality 

/ = TTf = Tff 
holds true for all of these three functions. I used the nomenclature: 

[Tf)y := f dx e-^f(x), 



Moreover, I define for integrable and bounded functions /, g the convolution of / and g 

as 

(f*9)(x):= / dy f{y)g{y-x). 
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(There are slightly differing definitions in the literature, hence I must clarify the 
notation.) The convolution theorem for Fourier transforms can be formulated by the 
equation 

P{f*g) = Pf ■ Tg. 
Application of this knowledge to the product a.\ f\ yields: 

F{anh)=Fan * Ph. (11) 
Integration by parts yields: 

where I use the notation z := A + iy and the definition 

F(z) := J dx e zx f'(x). 

(By the letter "/" I denote now the normalized Fermi-Dirac distribution; "/'" is its 
derivative.) The function F(z) is holomorphic on the stripe {Re e] — 1, 1[} within the 
complex plane. For purely imaginary arguments, its calculation is equivalent to the 
calculation of the Fourier transform of /'. This I let be performed by Mathematica 
and obtained the identity 



„. , — tnz 

F z) 



smh(—i7cz) 

at first only for purely imaginary arguments. However, the two functions are 
holomorphic on {Re g] — 1, 1[}, and so the equality holds on this larger set, too. 
One obtains: 

(// A ) y=- sinh ^ y _ iX ^y 

Using this information and the notation: 

7T Z 

H{z) := — —, — r (holomorphic on {Im e] — 1, If} , with H(0) = 1) 
sinh(7rz) 



the limit A — > of equation (11) turns into: 



{r«t)v = - iim / dx (J r ^ + yWx-iX) 

v ' 1 2tt a^o J A + ix 

One can integrate like 

//*oo 
dx function(x) = / (function(-x) + function(x)) ; 

after this, it makes sense to multiply numerator and denominator of the fractions 
with the complex conjugate of their denominators. Upon grouping the contributions 
according to their prefactors one obtains two summands. The first one reads after an 
integral transformation: 

(Fai) (-Ax + y)H(X(-x - i)) + {Fa{) (Ax + y)H(X(x - i)) 



oo 

dx- 



1 + x 2 
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— >■ 7r (J-ai) (y) (A —7- 0) with Lebesgue. 
The second integral is 



dx 



x 



\ 2 + x< 



{(J r ai) (—x + y)H(—x — iX) — (J-ai) (x + y)H(x — iX)} 



dx 



dx 



x 



X 2 + x 2 

X 

X 2 + x 2 



(Foci) {-x + y) {H{-x - iX) - H{x - iX)} + (12) 
H(x - iX) {(Fai) (-x + y) - (JFa{) (x + y)} . 



To the first line (12), the theorem of Lebesgue can be applied; because of the symmetry 
H(x) = H(—x), one obtains zero for it. For the treatment of the second line I use that 
Toli satisfies a Lipschitz condition of the form \function(a) — function(b)\ < const\a—b\. 
Hence, the convergence theorem can be applied also to the second line and one obtains: 

{Fai) (y + x)- (Tai) (y - x) 



dx- 



sinh(7rx) 

What remains to be shown is that this integral decays exponentially as a function 
of y. One can write J °° = + J xq and treat the two intervals separately. For example: 

{Fai) (y + x) 



dx 



■Co 



sinh(7nc) 



< const I dx \J-qli\ (x — (—y))e 



-ir\x\ 



Because of the exponential decay of Tati, this is bounded by the integral of the product 
of two exponential decays: One is centered in x = 0; the center of the other one is 
shifted by —y. This integral decays, as a function of y, exponentially. The integral 

(Fcti) (y - x) 



dx 



sinh(7rx) 

can be treated in the same way. 

What remains to be shown is the exponential decay of 
Xo {Fa l ){y + x)-{Fa l ){y-x) 
smh(7rx) 

as a function of y. To this end, one can write the absolute value of the integral as 

Xo {Fon) (y + x) - (Tax) (y - x) 2x 
dx 



2x 



< const ■L((Fa l )\ [y _ XOjy+Xo] ) 



sinh(7rx) 
: const ■ l(y) 



where for any function / defined on an interval D the notation 

I /(*)-/(») 



L(f) := sup 



x-y 



: x,y e D,x ^ y 



was used and where " f\ D ," is the restriction of a function to a smaller set. l(y) is finite 
for any value of y since 



(Foti) (z) = nc e 



-ix z -\z\cq 
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satisfies a Lipschitz condition as already noted earlier. (I finally used the Lorentzian 
shape of the function a/, Eq. (10)). For values of y with \y\ > Xq, one can use the mean 
value theorem in order to obtain an upper bound for it. The exponential decay of the 
derivative of Ton causes the exponential decay of l(y). 

Remark: I made use of the Lorentzian shape of ai only at a late stage. Therefore, 
I think that the proof works also for more general functions: One might take a function 
which is constant over a finite interval and convolute this with a Gaussian of total weight 
one and small width. This has a smoothening effect [5]. Taking linear combinations 
of such functions, one can approximate a quite arbitrary version of «/. At any rate, it 
would not be satisfying if the proof would only work for a Lorentzian. 



4. From the kernels to the current 

4-1- Representation of the current in terms of the kernels 
It was shown that for w < Wq: 

("K c " denotes the current kernel.) From the recursive definitions of the kernels j6], one 
can conclude on a general level (not only for the SIAM): 

• x a hermitian operator on the quantum dot space =>- K^'^x hermitian; 

• for any x: Tr (K^x) = 0; 

• for any hermitian x: Tr (^K^ n \x)^ real. 

For the SIAM, one can see by looking at the structure of the diagrams that any matrix 
of the form \a >< a\ ( with a — 0, t> I, 2) is mapped by every single order of the kernels 
to a linear combination of these four operators. Therefore, one can view the density 
matrix kernel and all of its orders as maps 

where V is defined as the four-dimensional real vector space spanned by the operators 
| a >< a\. Their image is always contained in the three-dimensional linear subspace 
U C V of the operators with vanishing trace. It is technically useful to introduce the 
operator 

L(w) := \k{w). 

The quantum master equation in the stationary limit can be written as 
L{w)p{w) = 0. 

The stationary reduced density matrix of the quantum dot, " p(w)" , is determined by 
this equation. 
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Since the image of L(w) has dimension three or less, the equation has a solution. 
L(0) is given exclusively by the second order diagrams [7] and, in the basis V\ : = |0 >< 
0|,f2 := | f><t 1)^3 := I i><4 |,^4 : = |2 >< 2|, it can be written as the matrix: 

/ ... a-(^io) a-(^io) \ 

a+(E w ) ... a~{E 2l ) 

a + (E w ) ... a~{E 21 ) 

\ a+(E 21 ) a + {E 2l ) ... / 

For this I used the abbreviations E\q = E, E 2 \ = E + U, where E and U are the 
energies appearing in the Hamiltonian of the isolated dot, Eq. Q, and the notation 
a ± ■= 

L(0)v i, L(0)t>2, L(0)v3 and are linearly independent. Only the assumption that 
a + (E w ) et cetera are strictly positive is needed for this. The linear maps V — > V are 
called "endomorphisms" ; the vector space of the endomorphisms is End(V). I want to 
use the map 

: End{V) ->■ End(V), 

<t>{E) (j2 x i v i] = X i E ( v i) + A ^4- 

\ i=l J i=l 

0(L(O)) is invertible, is smooth (it can be derived an arbitrary number of times). 
The set of the invertible endomorphisms, u {det ^ 0}", is topologically open in EndiV). 
Therefore, there is r > such that for all 

E G B r (L(0)) := {E' G End(V) : \E' - L(0)| < r} 

still 4>(E) G {det ^ 0}. ( All norms on EndiV) are equivalent, since it is a space of 
finite dimension.) This allows me to define 

F : B r (L(0)) V, 

F(E) := [inverse(4>(E)) EV4] — V4. 

For sufficiently small values of w, one finds L(w) G i? r (L(0)). Then, necessarily, 
the three operators L(w)vi,i = 1,2,3, are linearly independent and L(w)v4 is a linear 

combination of those three: 

3 

L{w)v4 = XjL{w)vj. 

8=1 

This implies that 

3 

F(L(w)) = x i v i ~ «4, L(tw) = 0. 

i=l 

The result obtained so far is that for, say, w < Wq, the space of the solutions to the 
quantum master equation in the stationary limit, 

L(w)p(w) = 0, 
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is spanned by F(L(w)), the dimension is one. The equation L(0)p(0) = is solvable by 
an operator p(0) with trace one. This implies that Tr {F(L(0))} ^ 0. The composition 
of maps 

Tr oF 

is smooth, and so one can conclude that there is r' < r such that 

for all E £ B r /(L(0)) : Tr(F(E)) ^ 0. 
As a consequence, the map 

/ : B r r(L(Q)) — >■ V, 

F(E) 

E i y 



Tr(F(E)) 

is still smooth. For w < w' Q , with fixed w' Q < w , one obtains that L{w) £ B r /(L(0)) and 

p(w) = f(L(w)). 
The equation for the current u I(w)" reads: 

I(w)=Tr K c (w)[f(L(w))). (13) 

4-2. Taylor series of I(w) in w = 

The above representation of the current in terms of the kernels shows that it is a smooth 
function of w. It has a well-defined Taylor series in w = 0; the coefficients of this 
power series are given by the derivatives of the current with respect to w in w = 0: 

Cfc = I/«( w = 0). 

(All odd derivatives vanish, since the current is also a smooth function of w 2 ; this is not 
relevant for the following arguments.) The Ck can alternatively be defined recursively 
without taking derivatives: 

c := I(w = 0), 

p (w) := c , 

(I-p k )(w) 
Cfc+i := am — , 

w— >0 W + 

Pk+i(w) := p k {w) + c k+1 w k+1 . 
This representation is equivalent to the definition by derivatives of I{w). 

4-3. Replacing the kernels 

Perturbation theory is applied [9] by cutting off the expansions of the kernels at some 
order. Upon defining 

2^( 2 ) i _L„„2n^(2n) 



K (c) (w) :=w 2 K$ + ...+w 2n K i 



(c) 
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the density matrix u p(w)" is determined by the equation 

K(w)p{w) = 0, 
and the current I(w) by 

I(w) = K c (w)p(w). 



In the same way as done for the exact current (13), this can be expressed in terms of 
the kernels as 

I(w) = Tr K c (w) [f(L(w))] for w < v%, 

where I defined L(w) := —%K(w). Although rather unusual within this text, I point out 
that the map "/" appearing here is precisely the same as the one which appears in Eq. 
@). 

For sufficiently small values of w, one can assume that both L(w) and L(w) are 
contained in B^(L(0)), and therefore: 

2 

/ [L(w)] = f [L(w) + w 2n b(w)} = f [L(w)} + w 2n B(w), 

where by u b(w), B(w)" I denote bounded functions. The reason is that / is smooth on 
B r '(L(0)) and thus the derivative of / is bounded on B^(L(0)). Moreover, one can use 

2 

the estimate 

|/(z) - f(y)\ <\x-y\ sup {\f(z)\ : z E B^(L(0))} for x,y E £^(L(0)), 

which is obtained from the mean value theorem. 

Analogously, the current kernel can be written as 

K c (w) = K c (w) + w 2n+2 b c (w) with bounded b c (w). 

Putting the pieces together, one obtains: 

I(w) = I(w) + w 2n+2 bj(w) with bounded bj(w). 

Inserting this expression for I(w) into the recursive definition of the coefficients of the 
Taylor series of I(w), one sees that the coefficients c , . . . , c 2n equal the corresponding 
coefficients of the Taylor series of I(w). Alternatively: 

l( 2k \ w = 0) = I^ 2k \w = 0) forA;<n. 

(The odd derivatives vanish, anyway.) 



5. Conclusion 



I showed a possible way to calculate the derivatives of the current across the SIAM with 
respect to the coupling parameter u w tl in w = 0. Implicitely, I showed the existence of 
these derivatives. Remark: In the same way as the kernels, the current is also analytic 
in w around w = 0. This can be seen by solving the quatum master equation in the 
stationary limit explicitely. A representation of the current directly in terms of the 
kernel elements is obtained from this. 
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6. Why a perturbation theory for the Anderson model? 

There is a very nice and straightforward interpretation of the current obtained from the 
second order kernels (e.g. Ref. [7j) in terms of energy conserving one-electron processes: 
If a, b are quantum dot states with particle numbers N(b) = N(a) + 1, and if the dot is 
in the state a, then electrons can tunnel from the leads to the dot in case their energy 
equals the difference Ef, — E a . On the other hand, unoccupied electron levels in the leads 
with this energy are needed for the inverse process. A difference between the chemical 
potentials of the two leads gives rise to an effective current " j( 2 )" of electrons across the 
quantum dot. Strictly speaking, the perturbation theory is the link between the object 
l( 2 \ which can be understood so nicely, and the final exact model of the current used 
within the real-time approach. This was my initial and main motivation for trying to 
see why exactly the statement of the perturbation theory is correct. The proof implies 
that indeed is closely related to the exact current. (On the other hand, I admit that 
the statement of the perturbation theory was perfectly plausible.) 

Recently, I noticed that it has been maintained that spin-flip scattering processes 
as assumed in Ref. [UJ are essentially contained in the SIAM with the conventional 
tunneling Hamiltonian, Ref. [3] (section 10.5.1). The coefficients of the spin-flip 
scattering processes obtained by the argumentation in this Ref. are of second order in 
the tunneling Hamiltonian. This might give rise to a zero bias anomaly in the transport 
which originates from the fourth order in the tunneling Hamiltonian. However, not 
any zero bias anomaly was seen in Ref. [H] by the fourth order theory, applied to the 
SIAM. The perturbation theory might be used for a serious attack on the statement 
that spin-flip scattering processes can be derived from the Anderson model. 
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